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ABSTRACT 


The solidified waste product produced from a fuel 
reprocessing plant has high heat generation' rates . It is 

therefore necessary to store these wastes for an interim 
period before final disposal . This is stored in well 
sealed containers to prevent the leakage of activity 
These containers are continuously cooled to dissipate the 
decay heat . 


For this purpose natural or forced circulation of 


air is used . 

For a 

given 

set of 

conditions 

heat 

generat 1 on 

rat 

e, heat 

transf er 

coefficient 

etc . 

have 

been 

estimated . 

The 

shape 

of storage canisters 

is 

cylindrical 


These calculations have been done for different orientations 
of canisters ,e.g. , vertical, jtiorizontal , vertical with sleeve 

C annular ) and cross flow across a tube bundle . 

After this the temperature distributions radially 
as well as along the height are evaluated using a finite 

difference scheme requiring a comprehensive computer 

programme . 

Finally accident analysis has been done i.e. 
i^ the entire system fails and air flow is completely 

choked ^ Then how long would it take for the temperatures 
to rise above critical limits , within which maintenance 
personnel' should be able to repair the system . 
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Chapter 1 


INTRODUCTION AND WASTE STORAGE FACILITY 

1.1 INTRODUCTION 

The epent fuel derived from a nuclear reactor ia 
reprocesaed to extract the unuaed fuel from it . The 
vaatea coming after reproceaaing are characterized by high 
heat generation ratea and are highly radioactive , 
eapecially if it ia of high active type . Therefore it 
becomea neceaaary to atore theae waatea for an interim 
period of time . During which a major portion of decay 

heat ia diaaipated ao that an inert and atable form ia 
available for the final diapoaal . Thia ia known aa 
conditioning of the waatea , after which it becomea eaaier 
to handle , tranaport and atore them aafely . After the 
interim atorage they are diapoaed to the ultimate diapoaal 
aitea with leaaer coata . 

The atorage facility where theae waatea are 

atored for an interim period haa to be deaigned very 
carefully . All the thermal , biological and atructural 

conditiona ahould be aatiafied by it . Thermally it ahould 
be able to. provide contlnuoua cooling of the product in 

order to- avoid melting , which may lead to an extremely 
dangeroua aituation . Biologically it ahould be able to 
contain radiationa coming out of waatea , and atructural ly 



£. 

it should be able to withstand tremors , etc. . To ensure 
the integrity o£ wastes and its shielding canister , 
constant surveillance is required . It is expected that 
thermal and radiation conditions are stabilized after a 


period of 

about 10 

years , 

but to be on 

saf e 

aide 

an 

additional 

period of 

about 

10 years is 

provided 

- 

Henc e 

the storage 

facility 

should 

have sufficient 

space 

to 

store 

the wastes 

generat ed 

over 

a period of 20 

years . 

[1. 

2] 


1.2 RADIOACTIVE WASTES [3.10] 

Basically there are following four kinds of 
radioactive wastes 

1. High active wastes ( HAU ) 

2. medium active wastes ( MAU ) 

3. Low active wastes ( LAU ) 

4. Transuranic wastes ( TRU ) 

HIGH ACTIVE UASTES 

HAU is a term used for wastes derived from 
partitioning fission products and transuranics from uranium 
and plutonium during the reprocessing of used nuclear 
fuel. It may be in the form of sludge , calcine , or other 
products £ienerated in treating liquid HAU . They 
require cooling and shielding both and great care is 
required while storing them .The radioactivity level of these 
wastes is 100 pci ml ^ and above. 
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MEDIUM ACTIVITY UASTES 

, intermediate 

MAU is the term used ^ 

r.^r,ulc6 cooling but 

level radioactive wastes which do not req 

* A but foC” which 

require shielding and can be safely store 

all Able and which 

disposal facilities at present are not av ' 

r Kaa+ venerating wastes 
do not fall within the category heat g 

^ rt ^ ^ a 


-5 

HAU . The radioactivity level of these wastes 


to 


100 ^ici ml 


LOU ACTIVITY UASTES 


LAU are defined as radioactive wastes 

^A are usually 

not require shielding and cooling ana 

. , radioactivity 

classified in terms of the concentration o 

they contain . Although LAU are those materials vi 

^^rtiiire treatment 

level of activity which do not often req 

And potentially 

before disposal , they may contain some low 

ti/Iaa Tbeir radioactivity 
hazardous cocentrat ions of radionucl iaea - J-n 

level is leas than 10 ^ pci ml ^ • 


TRANSURANIC UASTES 


Transuranic wastes are those elements 

isotopes beyond uranium in the periodic table 
produced almost entirely by the irradiation 
nuclear reactors . Characteristically these are ^ 
whose decay schemes follow the natural U and 


containing 
these are 
U 4 Th in 
emitters 
Th decay 
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schemes . Examples of these isotopes are 

plutonium^^^ , Armecium^^^ , Curium^^^ ’ , Neptuniun^^^ (very 
long half life , 2.14 million years ) . 

13 WASTES DURING REPROCESSING £3] 



Fissioning 

of nuclear fuel 

in 

a reactor 

creates 

highly 

radioactive 

materials 

which 

are 

almost 

all 

retained 

in the 

used fuel 

elements 

when 

removed 

from 

the 

reactor 


. Uhen the used fuel la reprocessed to recover unburnt 
uranium and plutonium , moat of the radioctivity is 

retained in the liquid wastes which are later converted 
into solids . 

Reprocessing of used fuel from a nuclear power 
facility of 1 GUe operating at a plant factor of 70 percent 

would generate 5000 to 7000 kg of solidified HAU . 

About 5000 liters of primary HAU solutions are 
used in the reprocessing of 1000 kg of used uranium 
fuel. The volume would be reduced to 1.1 liters by 
evaporation before storage in the stainless - steel' storage 
tanka . Prior to the actual solidification process , the 

volume would be further reduced by evaporation to some 

0.38 liters . 

Table 1.1 shows a typical inventory of radioactive 
wastes arising from reprocessing - 



TABLE 1.1 


Inventory of 

radioactive wastes during 

reprocessing 

Reprocessing wastes 

Waste material 

Plutonium 


generated 

cont ent 


( m^/Gg HM ) 

C Kg/Gg HM ) 

HAU 

500 

50 

LMAU 

1500 

15 

Combustible 

1845 

8 

Non - Combustible 

1500 

15 

Cladding spacers 

350 

15 

Dissolver residues 

50 

0 


In most of the fuel reprocessing plants , gaseous 
effluents generated are particularly important , even though 
some isotope decay quickly . In particular iodine -129 which 
has only a 0.8 percent fission yield has a half life of 
17 Myr , carbon -14 - 5500 years , krypton - 3950 days and 
tritium - 12 years . BUR fuel elements contain about one 
percent unburnt uranium -235 , 0.5 percent plutonium isotopes 
, mainly plutonium 239 , and the balance about 95 percent 


nonfissile uranium -238 . 
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1.4- ALTERNATIVES FOR WASTE STORAGE [ 4 , 5 ] 

Several deaigns have been considered for interim 
storage . Use of water ponds for providing interim storage 
for the solidified waste is an alternative . One option is 
to use air cooling on the basis of detailed techno - 
economic consideration , both air and water are 

acceptable .Uater of course is a better cooling medium 

because of its high heat capacity , as compared to that 

of air . However , consideration of corrosion , need for 
demineralised & chemically adjusted water & high make up 

water are some of the disadvantages . The other alternative 
uses natural convective air cooling for air cooled storage 

the waste canister is enclosed in a secondary container 
to ensure maximum safety against accidental release of 
activity in the event of canister failures . Cross flow of 
air with stack induced natural circulation is preferred 
However a back up system for forced supply might be 

required to cater for unforeseen eventualities leading to 
breakdown of natural air circulation . 

1.5 DESIGN OF SOLIDIFIED WASTE STORAGE FACILITY [l] 

The solidified radioactive wastes produced from 

reprocessing plant are stored in fully sealed , stainless 
steel canister . The air cooled storage utilizes decay heat 

& a suitably designed stack to provide the driving force 
required for the movement of air through the storage 
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vault . In effect the waste heat increases the air 
temperature causing an upward movement of air due to 

bouyant force . 

The design depends upon the rate at which 
containers containing solidified highly radioactive wastes 
are received from the waste immobilization plant . In fact 

the waste from the reprocessing plant is cooled for 900 

days before solidification . The liquid waste contains about 
675 ci lit.^ of ft , f activity generating about 2 watts 

lit.^ of decay heat . The specific activity of the 
vitrified product is around 13,500 ci lit.^ generating about 
40 watts lit.^. It is important to note that actinides have 
not been taken in to consideration in these estimates 

since they are relatively leas important within the time 
frame of interim storage . The actinides assume primary 
importance for long term storage and disposal condition , 
for instance in the Tarapur Storage Facility , which has 

been designed to store wastes generated over a period of 
20 years , actinides have not been taken in to account . 



To 

meet 

the design 

objectives 

a secondary 

containment 

called 

storage unit is 

provided 

to the storage 

canister 

. For canister simple 

cylindrical 

geometry 

is 

selected 

. Finned 

or annular geometry is not 

envisaged 

in 

view of 

low 

heat 

generation rate 

of the liquid waste 

due 

to moderate 

burn 

up of the spent 

fuel . For 

enclosing 

the 

canisters 

in 

the 

storage unit 

three arrangements 

are 


considered : 


-Axial stacking 





-Cortinatlon of two 

A onH third will reduce considerably 

Second 

overall storage positions - But it has got disadvantages of 

^ deficiency leading to uneven temperature 
poor heat transfer eru-'-* 

distribution in the csnistsrs . From heat transfsr point of 

view axial stackine 1= atoraga unit with 

Incrsaasd height . m aurv.illance facility 

at Tarapur two canlatera are atacked in a aingle storage 

, a. _ ftnlt is completely sealed 

unit . The storage unit 

the lid on . 

Thus in a storage unit C shown in 

there are three barriers for the outward 

activity - 

Float the borosllicate matrix in 

oxides are incorporated . 

rm w a-ie* hiah integrity of all welded stainless 
Second the ni«« 


by 


welding 


Fig. 1.1 ) 

movement of 


which 


waste 


steel canister & 

Third the ateel welded carbon steel 


secondary 


containment . 

For the purpose of surveillance activity the loss 
would msan partial or total loss of all ths thrss 

barriers . This would happen on the highly unlikely event 

I t u. , , , melting down & corroding both the 

of solidified waste mexLxa® 

primary i aecondary containment .If and when air borne 

activity ia detected m the cooling air , the storage 

nnit which has lost its Integrity can be Isolated . 

withdrawn , rescanned and restored to its position . 



323*65 nvn 
OO 




355* 6 mm 

b 


STORAGE CANISTER 


STORAGE 


Fie- 


1 .1 Schematic 


o£ storage un 
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1.6 NATURAL CONVECTIVE AIR COOLING STORAGE [13 

The air cooled storage is based on natural 
convection air cooling with induced draught . It utilizes a 


stack 

to provide 

the 

driving force 

for 

air 

movement 

through 

the storage 

vault 

. The system 

is 

self 

regulating 

& can 

compensat e 

for 

the 

changed load 

. As 

the 

decay heat 


reduces , both the air flow and air temperature decrease . 
Storage vau lt 

It is the building (shown in Fig. 1.2) in which 
storage units are arranged in a definite array , in order 
to optimize heat transfer efficiency . To facilitate this , 
storage vault is divided into two blocks and each is 

further subdivided into three compartments at Tarapur 
Storage Facility . It is necessary to isolate thermal 

effects from the load bearing structure and also isolate 
any water ingress . This is achieved by recourse to a 

double vault design . The inner storage vault is 

designed on thermal considerations & the design of 

external is based on structural and biological shielding 
considerations as well as protection from all credible 
accidents . 

The filling pattern of the storage vault may be 
progressive, or compartment wise . In progressive filling 
pattern • each compartment receives one storage unit 
successively . In this pattern the peak heat load in any 
compartment would not exceed about . 510 Kw , with an 



STACK 
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wise filling i.e. if a compartment is to receive 


One 


storage unit every alternate day 


Then the 


“axiffluu 


cumulative heat load would be about 820 Kw and 
load around 0.3 billion curies . Uhile progressive 


activity 


fill 


has the advantage of providing for lower heat 
activity loads in the vault structure , compartment 


ing 


and 


wise 


filling would provide isolation of radioactive 
in any one compartment . Each of these 


'entory 
filling 

pattern is further subdivided into two loadina « 

® patterns 

namely forward and backward , forward loadina 

* Bieana 

placing the unit in the direction of air flow 

backward implies the reverse . Backward has an 
that the fresh storage unit with maximum heat 
comes in contact with the cooler air while in 


loading the storage unit with maximum heat generation 
in contact with the hotter air . 


and 
advantage 
Seneration 
forward 

C 0111^3 


1.7 OUTLINE OF PRESENT WORK 


In this thesis convective heat 

titans fee 

coefficients have been estimated for different arrangem 

of storage units . The storage unit may be held VAt^+i 

^®*^tically 

or horizantally with free or forced convective 
cooling . It may also be held within a sleeve , in 
both the convection and radiation heat transfer 
simultaneously . 

In the steady state analysis of a storage unit 
the two dimensional temperature distribution has 

been 


air 

which 

Occur 
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obtained . The method used is finite difference . From the 

computer programme prepared the temperature variation can 
be obtained in a cicular cylindrical object , radially as 

well as heightwiae , subjected to forced or free convection . 

In the transient analysis part the time has also 

been regarded as a variable . It is very useful from 

accidental analysis view point , as it gives the results 
under extremely pessimistic conditions . Here , variation of 

atmospheric temperature , change in radial temperature 
distribution etc. with time has been obtained . 
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Chapt-er 2 

ESTIMATES OF CONVECTIVE AND RADIATIVE HEAT TRANSFER 

COEFFICIENTS 

For the estimation o£ convection heat transfer 

coefficient ( he ) different empirical relations have been 
used depending upon the type geometry and flow conditions 
. The various configurations are discussed one by one 
Throughout the medium is air . 

2.1 CALCULATIONS OF FLOW RATE AND PRESSURE LOSS THROUGH 
THE TUBE BANKS [6,9] 

Storage units are held vertically in a triangular 
pitch of 825 X 825wm -Hence the arrangement is staggered 
tube flow . ( Fig. 2.1 ) 

Given 

Air inlet temperature T = 25 C. 

1 

Maximum permissible air outlet temperature T _ = 85 C. 

00 2 

Total surface area for heat transfer is A = N.n.d.l 

where , 

N = Total number of storage units = 1848 units 

d = diameter of each unit = 35,56 cm. 

1 = length of each unit = 200 cm. 

hence 

A = 1848 X n X 35.56 x 200 
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Top view 


Storage uni 



Fig. 2.1 Arrangement of storage units in the storage vault 


1 

7 2 

= 4.13 X 10 Cm. 

2 

= 4130 meter 

From the heat balance equation 

Q ' “ x V ( 2 - 1 ) C2-1) 

where , 

Q = Total heat load = 4000 kW C However actual load 
ia 3060 kS'/, but for design purpose 4000 KW has been taken) 
C = Specific heat of air = 1.006 kJ kg 
m = Maas flow rate , kg s ^ 
hence 

1 _ 4000 

" 1.006 X ( 85 - 25 ) 

This gives 

ffl = 66.282 kg s ^ 

Density of air - 

p 

P - — ( From the gas law ) 

cts R T ^ 

oo 

where , 

P^= Density of the atmospheric air 
P = Pressure of the atmospheric air 
= 1.0132 X 10^ N m“^ 

R = Gas constant for air = 287 J kg ^ 

T = Temperature of the atmospheric air = 298 K. 

00 

hence 

p^= 1.0132 X 10^/( 287 X 298 ) 

= 1.185 kg m"^ 

Therefore volume flow rate of air is 

V = m/p 
00 

= 56.92 m^ s~^ 

= 2,01,324 m^ hr”^ 
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For pressure loss calculation The formula used is 


2 X f X G X I 

* "o rn3L3c t ^ _ V 

AP = — — X ( ) 

P ^ '^w b 


( 2 . 2 ) 


where , 

-2 -1 

G = mass velocity at the minimum area , kg m s 

max ' 

p - density evaluated at the free stream condition 

p, - average free stream viscosity , Pa. s 

D 

0 

^ = viscosity of air at the wall temperature (179 C.)» 

w ' 

f = friction factor 

= number of transverse rows = 66 
For the given arrangement 


f = 0.25 + 


0.118 

[(S^- d)/d3 


-0.16 


1.08 


X Re 


(2.3) 


where , 

Re = Reynolds number at the minimum flow area 
max 

The minimum flow area is 


A. =CS^-<DxlxN 
min n g 


For the 

given 

arrangement 


S 

n 

825 mm . 



d 

Diameter of the storage 

1 

1 ength 

of the unit 

= 200 

N 

8 

number 

of gaps for 

the 

for the 

given 

arrangement 

= 32 

hence 




A . = 

min 

C82.5 - 

35.563 X 200 x 

32 X 

= 

30.0 m^ 




.-4 2 



IS 


And inaaa velocity at the minimum area 


G = m/A , 
max min 

= 2-21 Kg/m^aec 

flow velocity at the minimum area 


V = V/A , 
max min 


= 56.92/30.0 


= 1.897 m a 


-1 


Also 


P 

= 1.185 

CO 

-5 


= 2 X 10 Pa. a 


-5 


= 2.4 X 10 Pa. 


w 

hence 

Re 


(at 25 C. ) 
(at 179 °C. ) 


= Cp X V X d 

CO max 1 


max 

= Cl. 185 X 1.897 X 0.3556 3 / 2 x 10 
= 39968 

And hence friction factor is 


-5 


f = 


0.25 + 


0.118 


£C 82. 5-35. 563 /35.56] 


1.08 


X (39968) 


-0.16 


= 0.062 
Therefore 


, C 2 X 0.062 X 2.21 x 663 , 2.4 .. 

Ap = xC 3 


0.4 


= 15.66 N m 


1.185 

-2 


2.0 


-2 


= 1.596 kg m 
= 1.50 meter of air 

Thia pressure loss is very small in comparison 
head developed by the stack , which is around 


2.4) 


to 


100 meters 
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o£ air hence it does not affect the stack height much . 

2.2 CONVECTION HEAT TRANSFER COEFFICIENT FOR FLOW ACROSS 
STORAGE UNITS C 6 ] 

The formula to be used for the estimation of 
average heat transfer coefficient for flow across the tube 
bundles is (shown in Fig. 2.1) 



Uhere , 

-2 -1 

h = Convection heat transfer coefficient , U m K 
d = Tube diameter , m. 

= Thermal conductivity of air at film temperature ,U m 
U = Free stream velocity of air , m a ^ 

Pr = Prandtl number at the film temperature 

-3 

= Density at the film temperature , kg m 
/-<£ = Dynamic viscosity at the film temperature , Pa. s 
c and n are constants for a particular arrangement . 
For the given arrangement 
C = 0.523 





n =0.56 

The film temperature is 
T^=CT + T )/2 

£ W 00 


where 


h 


T 


w 

T 


00 

ence 


Temperature of 
Temperature of 


storage unit wall = 179 *^C. 

T. + T ^ 

. , ... in out 

the surrounding air = 15 


T = (25+85)/2 = 55 °C. 

<£} 

T. = 117 °C. = 3SG K 

r 

At this temperature 

= 2.4 X 10 ^ kg m ^ 
= 0.03 U m"^K”^ 

Pr = 0.80 


and 


Re = 39968 
max 

hence 


h 35.56 1 

100 0.03 


same as in section 2.1 ) 

0.56 1/3 

0.523 X C 399683 x C 0 . 83 


this gives 

h = 15.46 U/m? °C. 

This is the convection heat transfer coefficient 
for flow across the storage units . 


2.3 FREE CONVECTION HEAT TRANSFER FROM HORIZONTAL AND 
VERTICAL CYLINDERS [6.7] 


Horizontal cylinder 


The correlation for free 


convection 


heat 


transfer 



from a horizontal cylinder is (shown in Fig. 2.2) 

-jl/6 

Gr Pr 


= 0.60 + 0.387 


[1+ C0.559/Pr)^^^^3^^^^ 


( 2 . 6 ) 


for 10 Gr Pr < 10^^ 

where , 

Nu = Nuaselt number 
Pr = Prandtl number 
Gr = Graahof number 

All properties are evaluated at the film 
temperature . 

Given 

diameter of cylinder d = 0.36 meter 

length of cylinder 1 =2.0 meter 

Ambient temperature of air = 30*^C. 

Rate of heat generation in the cylinder q’’' = 12.96 K.H'm^ 


1 =2.0 meter 


T = 30 C. 

CO 


here h is calculated by trial and error method 


assume 

h = 10 

U/mTK 

, he 

q’ ’ ’ X .,■ X 

4 

d^ 1 = 

h X n 

d 1 : 

where all 

the terms 

are 

this gives 

T = 

T + 

q . I 1 

w 

CO 

4 h 

Putting numerical 

values 


W CO 


(2.7) 


T = 145.2 C. 
w 

therefore 

T. = (T + T )/2 = 87. 6°C. = 360.6 K 
f ' w co" 

At this temperature 

^ - jKiV)emdtrc viscosity of air = 21.85 x 10 m a 

= 0.03 U m"^ K"^ 
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8 . . 

Pr = 0.691 and Gr Pr = 2.04 x 10 C which is in the 

given range ) 

Putting these values in the above correlation 
Nu = 70.17 
or 

h d/K = 70.17 

putting d and K into this it gives 
h = 6.12 U 

This value is much different than the assumed one hence 

repeating above calculations with this h 

Now 

= 124.13° C. = 397.13 K 
Pr = 0.686 

-6 2 -1 

i> = 25.92 X 10 m s 
K = 0.034 U m”^K“^ 

Gr Pr= 2.14 x 10° 


and 

putting 

these values 

into 

above 

correlation we get 


Nu 

= 71.14 






This 

gives 






h 

= 6.74 

U mT^K"^ 





Uhich is 

quite near 

to the 

last 

value , hence 

the 

value 

of 

convection 

heat 

transfer coefficient 

for 


— 2 -1 

horizontal cylinder is 6.74 Urn K 

Vertical • cyl inder 

For vertical cylinder with same dimension as for 
horizontal cylinder the correlation used is (shown in Fig. 
2.3) 




Fig. 2.2 Free 

cyl ind 



Fig. 2.3 F 



Z-Om 
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Gr Nu = 

X X 


C q”x 
“2 


= Gr 


K V 


C2.8) 


For laminar range : 
Nu 

a.nd 


* 1/5 

= 0.60 X CGr X Pr) ' 

X X 


for 


h = 5/4x h T 
av. x=l 

?or turbulent range 
Nu 


= 0.17 X (Gr*x Pr)^^^ 

X 


for 


ind 


5 * 11 

10^ < Gr < 10"" 

X 


13 ^ 16 

2x10 < Gr < 10 

X 


h 

av • 


X=1 



rhere , 





Gr 

X 


Grashof 

number at x 


Nu 

X 

= 

Nusselt 

number at x 


* 

Gr 

X 


Modified 

Grashof number at x 

ft 

= 

1/T^ 



K 

= 

Thermal 

conduct ivity 

of the air 

M 

= 

Dynamic 

viscosity of 

air 

X 


Distance 




ere again trial and error method will be used 

2 

33ume h = 10 U/m.K 
id 


^ d^l = hxfTdlx(T - T) 
4 ^ V oa' 


= 87.6 C. = 360.6 K 


givea 

( as was done for 

rlinder ) ■ 

: this temperature properties of air are 

V = 21.85 X lO"^ m’^s"^ 


horizontal 


K =0.03 U m“^K~^ 
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Pr = 0.691 


and 


ft = 1/T^ = 2.77 X 10“^K ^ 

Putting these values in above correlations 

* 12 4 

Gr = 1.84 X 10 X 

X 

* 11 

Now up to Gr = 10 


i . e . up to 


X = ( 10^/ 1.84x10^^) 


the flow is laminar 
1/4 

= 0.48 m . 


= 48.0 cm. 


and 

after 

this the 

flow 

is turbulent 

Now 

with 

* 

this Gr 





X 



For 

laminar range 

h 

X 

A , 0-2 

= 4 . 9 / X 

and 

f or 

turbulent 

range h 

6.14 U/m?K 

Now 

for 

the whole 

range 

the average 


U/m?K 


is 


av . 


0.48 
0 ^ 


aminar 


1.0 

dx + J' 

0.48 


^x turbulent 


= 3.40 + 3.20 
= 6.60 U mT^K"^ 


This value of 


IS 


quite 


di f f erent 


from 


assumed one hence repeating above calculations with this 
= 117. 3°C. = 390.3 K 


= 25.45 X 10 ^ m^s ^ 
K = 0.034 U m“^K"^ 


Pr = 0.686 


* 17 4 

Gr = 1.15 X 10 X 

X 

Now flow is laminar up to x = 54.0 Cm. 

0.2 


for laminar range 


h^= 4.83/x 


the 

h 
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& for turbulent range 5.38 

“2 -1 

on averaging above two h = 6,20 U m. K 

Again repeating above calculations with new h we get 

h = 6.17 U 
av . 

This value is quite close to the previous one hence 
for vertical cylinder the convection heat transfer 
coefficient is 6.17 U m.^K ^ 

2.4- FREE CONVECTION HEAT TRANSFER THROUGH A VERTICAL 
ANNULAR SPACE 


The empirical correlation used for this case (shown in 
Fig. 2.4) is 


Nu^ = 0.42.(Gr^.Pr)^'^^.Pr°'°^^.(l/<5)~°-^ 


where , 


(2.9) 


<5 = annular gap = 5 Cm. 

and all other terms have usual meaning , all the 

properties are evaluated at the film temperature . Here 

again trial and error method will be used 

“2 —1 

assume h = 10 U m. K , this gives ( from the heat 

balance equation used in section 2.3 ) 

T = 145. 2°C. and 
w 

T^= 87. 6°C.= 360.6 K 


The Grashof number is given by 




g^(T - T )^‘ 


1 > 


( 2 . 10 ) 




where all the terms 


are 


same as before 


Evaluating properties of the air at and putting 

above we get 
Gr = 819612.1 

and 

Gr-.Pr= 819612.1 x 0.691 = 566352 

o 

Putting these values in the above correlation 
Nu^= 3.80 = h<5/K 

O' 

hence 

h = 2.35 U 

Now with this h repeating above calculations we get 
= 275. 1°C. = 548.1 K 
Gr = 553322.1 
pr = 0.676 

and 

Nu^= 3.42 
which gives 
h = 3.0 U 

Again doing the same thing with this new h 
Gr^= 668171.6 

o» 

Pr = 0.683 
Nu^= 3.59 
this gives 

h = 2.94 U 

This value is quite near to the previous one hence 

— 2 -1 

for vertical annular space h = 2.94 U tn K 




into 



2.5 RADIATION HEAT TRANSFER THROUGH A VERTICAL ANNULAR 
SPACE [ 6 ] 


An alternative arrangement in the storage vault 
may be such that the storage units are held in a 
vertical annular apace , and there is a cross flow of air 
through the banka of storage units . After this the air 

is made to pass through annular space of individual 


units . In this 

way 

there 

will be two 

► convection 

heat 

transfers - 







- First due 

to 

cross 

flow of air 

through 

the 

storage 

units, it occurs 

at 

the 

outer surface 

of the 

sleeve 

- 

- Second due 

to 

free 

convection of 

air 

through 

the 


annular space between the storage unit and its sleeve . 

In addition to these there is a radiation heat 

transfer between the sleeve , storage unit and surrounding . 

Hence the total heat transfer coefficient ( h , ) in this 

^ total 

case is 


^total 
where , 


h 


c 


+ h 

r 


( 2 . 11 ) 


h = Free convection heat transfer coefficient for the 
c 

annular apace 

h^ = Radiation heat transfer coefficient for the storage 
unit 


h^ has already been found in the previous 


section 

. In order to 

get 

h , 
r ’ 

radiation 

analysis 

of 

the 

annular 

space (shown 

in 

Fig. 

2.5) is 

carried 

out 

The 


case is treated as, one surface held inside 


an 


enclosure 








Fig. 2.6 Radiation network 


where 


= net heat loaa from aurface 1 

^2 ~ net heat loaa from aurface 2 

= net heat loaa from aurface 3 

= abaolute temperature of aurface 1 

= abaolute temperature of aurface 2 

= abaolute temperature of aurface 3 = 303 K 

^3 emmiasivltiea of three surfaces 

Fj^j= shape factor from aurface i to aurface J 

_8 _2 -4 

O' = Stefan - Boltzmann constant = 5.669 x 10 U m K 


Shape factors and F 22 given dimensions 
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Fig. 2-7 Final network 
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in this network 

1-e. 




= 0.295 


1 12 


= 0.440 


= 0.231 


now from the radiation network 


qi + Qa = 0 


4 4 

- (tJ - tI ) 


( 2 . 12 ) 


The energy balance for the cylinder and sleeve 


-^gen = '^llrad 


(2 . 13) 


- ^2^ " ‘^llrad = Mouf'^2-^^2 " ^3^ 


(2.14) 


where 


*£ = film temperature = 


T +T 
1 2 


h . = convection heat transfer coefficient for the 
c ' in 

surface 1 and the inside surface of sleeve 


= 2.94 U m”^K“^ 


( calculated in section 2.4 


h I .= convection heat transfer coefficient for the 
c ' out 

outside surface of sleeve 


= 15.64 U m“^K“^ 


(calculated in section 2.2) 


q^en “ heat generation rate in a storage unit 
= 2164.5 U 


qil d7 heat loss from surface 1 due to radiati 
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= ambient 
on adding equns . 
of 


temperature = 30*^C = 303 K 

2.13 and 2.14 and putting 


the value 


q + h I . .A,. (T,-T_)/2 = 

gen c ' in 21 2 




putting numerical values into it 

2164.5 + 4.25 (Tj^-T^) = 3.32 (T^- T^) + 45.2 (T2-303) 
on simplifying this 

15860+0.93 

^2 “ 46 . 13 


C2.15) 


now put 

the values 

of , 

gen 

*^1 1 rad 

and 

1 2 into eqn . 2.13 

r t"* 

<y X [ - 

0.0113 - 

287.1 - 

3268345 

^1 - 

1.39 X 10^° 1/0.966 

+ 3.25 

- 1141.45 = 

2164.5 





“* 8 

putting <>• = 5.669 x 10 and simplifying above equation - 

5.87 X 10“®tJ - 6.63 x 10 "^°tJ - 1.685 x 10“^T^ + 3.06 = 4121.7 

— C2.16) 

Solving above equation numerically for 

= 463.4 K , and from eqn. 2.15 
= 353.15 K 

‘^llrad= 1^32.40 U . 

^ 1732.40 

^r SpiT-^-T^ 


hence 


and 



This is the radiation heat transfer coefficient 

for an annular space , and its value is around two times 
of that due to convection heat transfer . The value of h 

r 

may vary with the variation of annular gap and 
emmissivit i es of the surfaces . These variations are shown 
in Tables 2.1 and 2.2 


TABLE 2.1 


Variation of heat transfer coefficients with annular gap 


ANNULAR 

GAP 

C<5 ) , cm 


CONVECTION HEAT RADIATION HEAT 

TRANSFER COEF. TRANSFER COEF. 






TOTAL HEAT 
TRANSFER COEF. 

^^total\ura“^K^ 
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TABLE 2.2 


Variation of radiation heat 
emmissivities of surfaces ( 


trans f er 
and ^ 


coefficient with 
for a constant <5 



2.6 RADIATION HEAT TRANSFER IN FLOW ACROSS STORAGE UNITS 

In this case radiative heat transfer occurs 
between the storage unit and the surrounding concrete 
structure .The surface of neighbouring storage unit does 

not affect the radiation heat transfer , because they are 

at the same temperature . The units are held in a 
triangular pitch of 825 x 825 mm. (shown in Fig. 2.1). This 

case is treated as , one surface completely enclosed inside 

the other surface . 


The network of radiation heat transfer is (shown in Fig. 

2.8 ) 
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wher 6 , 

" temperature of the cylinder aurface 
= 452 K ( from section 2.2 ) 

T 2 = 323 K ( assume ) 



= emmissivity 

of surface 1 

= 0.6 

"=2 

= emmissivity 

of surface 2 

= 0.6 


= area of surface 1 

= Ti . d . L 

d 

= diameter of 

the cylinder 

= 0.36 m 

1 

= length of 

the cylinder 

= 2 . 0 m 


= area of surface 2 


= 2.2 [ 


0.825 X 0.825 y3 
2x2 


± .4.0.36 ] 



now for the network of Fig. 2.8 



SSa 


4 4 

o- ( Ti - ) 


1-e 


4^1 


1 ^ 1 


1-G, 


(2.17) 


^ 2^21 ““ 2^2 


and the , shape factor ^21 " » since surface 2 sees 1 only 

hence equn . 2.17 reduces to 

4 4 

<y.C ).A^ 


2 2 


putting numerical values Into it 

= 872.0 U 

hence , 

. 872 

= 3.0 U 

Thl« la the radiation heat tranafer coefficient , while the 

-2 - 1 

h for thla caae uaa 15.46 U m K , hence the total heat 

c 

“2 “1 

tranafer coefficient in thla caae ia 18.46 U m K . 


2.6(a; RADIATION HEAT TRANSFER WITH PERFECTLY BLACK SURFACES 

From section 2.5 it is clear that in caae of 
annular space the radiation heat tranafer coefficient rises 
rapidly with the increase of emmiaaivlty . Therefore it la 
possible to achieve a higher heat transfer rate by making 
the surfaces perfectly black. Under these conditions radiation 
analysis of both the previous cases ia discussed below 
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ANNULS SPACE 


In this caae (shown In Fig. 2.5) all the parameters 
leaving emmlsslvltl es are same as In section 2.5 . The 
emmlsslvltl es In this case are 


= 1.0 ; 


, - 1.0 ; 


.= 1.0 


with these values various resistances of the radiation 
network (shown In Fig. 2. 7) become 


1 - 


^lAi 


0.0 ; 


^ 1^12 


4 m a 
2 2 


= 0.0 


hence net radiation heat transfer is 


wher 6 , 


qi= ^ ( tJ - ) 


A.= Area of surface 1 = 2.26 m^ 


(2.18) 


Using same procedure as in section 2.5 to 


evaluate and , we get 


-9 3 

k ' nr* 


1.29 X 10 'T^ - 9.2 10 'T^ - 3.7 x IO'^tJ + 2.83 = 5096.84 


15860 + 0.93T 


-(2.19) 


( 2 . 20 ) 


solving equn. 2.19 numerically and putting the value of 
thus obtained into equn. 2.20 , we get 
T^= 417.6 K = 144.6 °C 
1 ^= 352.2 K = 79.2 °C 
now from equn. 2.18 


q, I 1935.8 U 

1 ' rad 


therefore radiation heat transfer coefficient is 
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1 ' rad 


-K--ra ' ^3 


T*= Ambient temperature = 30 C 


= 7.47 U m"^K ^ 

This la around 250?: of convection heat tranafer 

coefficient . Thia ahowa that radiation contrlbutea a major 
portion of total heat tranafer in thia caae . 

FLOU ACROSS STORAGE UNITS 
In thia caae 

1.0 ; 1.0 and all other parameters are same 

as in section 2.6 . Various resistances of the network 

of radiation heat tranafer (shown in Fig. 2.8) are 




= 0.0 


^ 2^2 


= 0.0 


^ 2^21 


1 

a: 


and hence 


4 4 

q^= O'. A2 . C - T^ ) 


( 2 . 21 ) 


putting same numerical values in this aa in section 2.6 

q = 5.669 X lo"® x 0.975 ( 452^ - 323^ ) 

^1 


1705.5 U 


hence 


1705.5 
r= . (T^-'rp- 




= 5.85 U m’’^K~^ 

Thia la only about 
tranafer coefficient in thia caae . 


of 


convection 


heat 



3 


Chapter 3 


STEADY ST ATE TEMPERATURE DISTRIBUTION 


Here ateady atate temperature distribution 
a solid circular cylinder , generating heat at a 
rate will be eatimated . The method used is 
Difference Method . 


inside 

constant 

Finite 


3.1 BASIC EQUATION [6] 


The basic heat conduction equation in 
coordinates is : 

1 S T . 








<>r 


r ^ r 


r^ S 0^ 


S 2 


q 

T 


1 » T 

Ot <9 T 


cylindrical 


(3.1) 


Uhere , 

T = temperature 
r = radius 
0 = base anfile 
z = height 
T = time 

q* ' * =s rate of heat generation 
K « thermal conductivity 
ot « thermal diffusivity 


For ateady state ^ T 

S T 


0 
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Chapter 3 


STEADY STATE TEMPERATURE DISTRIBUTION 


Here steady state temperature distribution 
a solid circular cylinder , generating heat at a 
rate will be estimated . The method used is 
Difference Method . 


inside 

constant 

Finite 


3.1 BASIC EQUATION [6] 


coordinates 
.2 


^ T 


The basic 
is : 

1 ^ T .3 


heat conduction 


+ — -= 




r S r 


r^ S 0^ 


a z' 


equation in 



1 ^ 

cn a r 


cylindrical 

(3.1) 


Uhere , 

T = temperature 
r = radius 
0 = base angle 
z = height 
T = time 

q’’’ = rate of heat generation 

K = thermal conductivity 
ot = thermal diffusivity 


For steady state a T 

a T 


0 
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3.2 , FINITE DIFFERENCE METHOD [11,12,13.14] 

In thia approach , the problem domain ia 
' diacretized ' ao that the dependent variablea are conaidered 
to exiat only at diacrete pointa . Derivativea are 
approximated by differences resulting in an algebraic 
representation of the partial differential equation . The 
nature of resulting system of algebraic equations depends 

upon the character of the problem poaed by the original 

partial differential equation . 

Here network lines are given by 
r^ = (i-1) X Ar and 

= Cj-1) X Az i , j = 1 , 2 , 3 

Quantities Ar and Az are increments in r and z 
directions respectively . The points of intersection of the 
network are called nodes . The network and nodes are shown 


in Fig. 3.1 



1 


z 


Figure - 3-1 : Deacret izat ion o£ the problem 




K \tf 
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various difference equations are : 
Forward difference 
T 


S T 


S r 


T. . . 

1 + 1 , j 


T. . 
1 . J 


i + 1, j 


Ar 


Backward difference 

^ T ^1 

^ r r i , j 


T. . - T. , . 

> j i 1 , j 

Ar 


Central di £f er ence 


^ T ^ , 

= T 

S r r 


1 » J 


T - T 

2 Ar 


Out o£ the above three di££erencee the central 
£ormula will be uaed . 


di££erence 


- T 1 

a ~ "''’li.j 


T. . . - 2 X T. . + T. . . 

i-1 . J 1 . J 1+1 . J 

Ar^ 


similarly 
T 

2 ^ 




a z 


. . . - 2 X T, . + T. 

i , J-1 i-, J J- > J + 1 

Az^ 


= T 


zz 


i > J 


Putting these expressions in equation 3-2 and simpli£ying 
2 ' “ ■^x.J ‘ ^2 * " ' ■ 


^ " ^i+1, j ^ ^ TcI 


i-1) 


) - h C 


(3.3) 


Where , C = -ijj- ^ h = aJT , K = AZ. 

equation 3.3 which involves £ive node points , is 
evaluated at each o£ the node points of network formed 
in this way n linear algebraic equations are formed if 
there are n number of node points in the network . In 
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For boundary no. 1 
S T 


S r 


= 0 


1 . e . 


( at r = 0 ) 

. 

^ » J 2 , j 


For boundary no. 2 


( at z = 0 ) 


or 


d T 

- K X — - = h X ( T. - T . ) 

^ z 1,1 amb . 


T - T 
12 ■^11 

- K X 12±- il_L_ = h X ( T, . - T . ) 

Az i , 1 amb . 


For boundary no. 3 
S T 


( at r = R ) 


- K « ^ = h x C ) 


or 


- K X 


T - T 

n-H,j n,j 

Ar 


= hx(T^, .-T .) 

n+1 , j amb . 


For boundary no. 4 


5 T 

- K X = h X C T. . - T ^ 

6 z i,n+l abm. 


( at z = H ) 

) 


or 


- K X 


T - T 

i,n+l i,n 

Az 


^ ^ ^ ^ ^n+1, j “ ^abm. ^ 


C3.4) 


C3.5) 


(3.6) 


(3.7) 


where , 

T , = Ambient temperature 

amb . 

In this way there will be n linear equations 
involving n unknown temperatures , these will be solved to 
obtain the temperatures . The methods of solutions are 


described below . 
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3.3 METHODS OF SOLUTION [ 11 , 14 ] 

There are numerous methods of solving n linear 

equations in n unknowns . These can be written in matrix 

form as 

A . X = B (3.8) 

Uhere / matrix A is the coefficient matrix 

matrix X is the column matrix of unknown 
temperatures , and matrix B is the column matrix 

whose elements are the right hand members of linear 

equations . 

Ilethods of solutions can be classified as direct 

or iterative . Direct methods are those that provide the 

solution in a finite and predetermined number of 

operations using an algorithm that is usually complicated , 
Iterative methods consist of repeated application of an 
algorithm that is usually relatively simple . They yield an 
answer only as a limit of a sequence . The number of 

iterations required to obtain solutions that no longer 

change to a specified number of decimal places usually 
cannot be determined in advance . 

3.3.1 DIRECT METHODS 

cramer * s rule 

This is one of the moat elementary methods 
introduced early in the study of algebra . Unfortunately , 
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the number o£ operations required in the algorithm is 
approximately (n+1)! , for n unknowns . For more than about 
three equations the use of cramer’s rule becomes 
impractical and is not recommended . 

Gaussian elimination 



It 

is a very 

useful and 

efficient tool for 

solving 

many 

aysteih of 

algebraic 

equations . Although it is 

one of 

the 

earliest 

methods 

proposed for solving 


simultaneous linear equations , it remains among the moat 

important algorithms in use today . 

The solution vector remains unchanged if the 

equations are multiplied or divided by a constant and if 
any equation is replaced by the sum or difference of 

that equation and any other equation . These facts are 

utilized in the elimination strategy . 

In Gaussian elimination , the objective is to 
transform the system into an upper triangular array by 

eliminating unknowns through algebraic operations . All the 

non zero coefficients will then be on or above the main 
diagonal . At that point the n th equation has only one 

unknown; the (n-1) th equation , only two unknowns; etc 
Back substitution , starting with the n th equation , readily 
provides the solution . The technique consists of two 

phases : the elimination of the unknowns and back 
substitution to obtain the solution . Loss of accuracy can 
occur in the algorithm due to round off errors . The use 



4 ? 


of partial pivoting to avoid division by zero and to 
improve accuracy will be included as part of the 
algorithm . If the elements of the coefficient matrix vary 


greatly 

in size 

, it 

is 

likely 

that accuracy 

can 

be 

further 

improved 

by 

the 

use of 

scaling . 




El iminat i on Method f or Tridigonal Matrices 

Gaussian elimination is a general procedure that , 
in principle , will provide a solution , if one exist , for 
any system of linear algebraic equations . For systems 
whose coefficient matrix contains a specific pattern of 
zeroes , the elimination scheme can be simplified to avoid 
pointless operations on the zero elements . Such specialized 
system of algebraic equations occur very frequently in the 
numerical solutions of partial differential equation . The 
moat common of these are tridiagonal and pentadiagonal , 
corresponding to equations that contain three or five 
unknowns respectively . 

Other direct methods 


Direct methods 
algebraic equations that 
Guassian. elimination do 
completely general . That 
the algebraic equations 
difference equations and 


for solving certain system 
are significantly faster 
exist . But none of them 
is , they are applicable only 
arising from a special class 
associated boundary conditions 


of 

than 

is 

to 

of 
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3.3.2 ITERATIVE METHODS [ 14 ] 

Iterative methods are particularly useful for 

system in which roundoff error may be a problem in that 

such schemes , when convergent , can be continued until 

changes in the solution have been reduced to some 

prescribed tolerance . Furthermore , iterative procedures can 

easily take advantage of the sparse nature of the 

coefficient matrix that are associated with many of the 

physical problems .( A matrix is sparse when a high 

percentage of the entries are zero . ) On the other hand , 
iterative procedures are certain to converge only for 
systems having "diagonal dominance” . Fortunately , many sets 
of linear algebraic equations originating from physical 

systems exhibit this diagonal dominance . 

Iterative methods can be further broken down into 

point C or explicit ) iterative methods and block ( or 
implicit ) iterative methods . For point iterative methods , 
the same explicit algorithm is applied iteratively to 
compute the unknown functions , whereas in block iterative 

methods , subgroups of unknowns are singled out for the 

solution by elimination schemes in an overall iterative 
procedure . The most common applications of the block 
iterative method are special classes of equations arising 
from the numerical solutions of elliptic partial 
differential equations . 
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Gauaa Si ed&l iteration 

It is perhaps the simplest of the commonly used 
iterative methods . Uhen it can be used , the procedure for 
a general system of equations is to 


1) 

make 

initial 

guesses 

for all 

unknowns ; 


2) 

solve 

each 

equation 

for 

the unknown 

whos e 


coefficient is largest in magnitude , using guessed values 
initially and the most recently computed values thereafter 
for the other unknowns in each equation ; 

3) repeat the solution of the equations in this 

manner until changes in the unknowns become smaller than 

prescribed tolerance . 

Actually a careful check on step 1 reveals that 
a guess for the unknowns having the largest 


coef f ieient 

in the 

f irst 

equation is 

not 

required 

. It 

should be 

noted that in 

carrying out 

step 

2 , it 

is 

necessary 

to determine a 

di f f er ent 

unknown 

from 

each 

equation . 

In most 

instances , when 

the 

Gauss 

si edel 

iteration 

converges 

, a different unknown can 

be found 

to 


have 

the coefficient 

largest 

in 

magnitude 

in 

each equation 

. Uhen 

this does not 

occur , 

step 

2 should 

be 

modified to 

permit 

all the unknowns to 

be 

determined 

from 

1 the system 


of equations . 

A more formal statement of the algorithm can be 
given if the equations are first ordered , if possible , so 
that the coefficient largest in each row is on the main 
diagonal . Then , if the diagonal elements are all nonzero ( 
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it must 

be 

possible 

to avoid 

a zero 

on the main 

diagonal 

if 

a unique 

solution 

exists ) 

, and initial 

guesses 

for 

the unknowns have 

been made , 




where it is understood that the moat recently computed 
values of x are used on the right hand aide ; that is , 

the variables are continuously updated . Uhen old values of 
X are used on the right hand side until all n variables 

have been computed , the iteration procedure is known as 
Jacobi itoratlori .In the Jacobi procedure , the variables 
are finally updated before the algorithm is repeated for 

another iteration . The Gauss Siedel algorithm converges 

more readily and should always be used . 

3.4- TEMPERATURE DISTRIBUTION 

The method used here is Gauss - Siedel iteration 

method for solving the linear algebraic equations . 

Following are the dimensions of the storage unit 

Length 1 = 2.0 meter 
Diameter d = 0.36 meter 

in fact these are the dimensions of the storage unit 
of Tarapur Waste storage Facility , so that results 

obtained can be verified later . 

other parameters used in the computer programme are 

= 10 U mT^K"^ 


Convection heat transfer coefficient 


h 



Rate of heat generation in the cylinder q’’’ = 12.96 kU m 

Thermal conductivity of the material 


K = 1 U m“^K"^ 


Ambient temperature 
Length of cylinder 


T = 60 C. 

CO 


1 = 2.0m 

Diameter of cylinder d = 0.36 m 

In fact the value of convection heat transfer 
coefficient may be changed in the computer programme in 
this part calculations have been done using four values 
of h in order to get a vide range of results . 

For the verification of temperature distribution 
obtained by finite difference method following equations 
are used : 

dq = hxdA>:(T - -T.. ^ ) 

surface ambient 

integrating 


q = Jh X C T 


. -T .. ,)xdA=q’’’x volume 

surface ambient 


this gives the total heat generated within the cylinder . 
Uhere , 


Hence 


q = Total heat generated in the cylinder 


oRiS.RY 







R is the cylinder radius 
dz is elemental height 

H 

q=Jhx2xnxRx ( ) x dz 

° (3.9) 

To evaluate above integral over the whole surface 


of storage unit 


1 y3 Rul 6 is 


used 


In 


this 
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method. Z is discretized into finite number of increments , 
each of size Az such that H = n x Az . The rule can be 
stated as follows : 

zl+n . Az 

J fw dz = ^ X [ + 4/2 + 2/3 4 /, * ... * 4 /^ 4 ] 

zl 

(3.10) 

Uhere /^ = / ( zl + (n-1) Az) 

From above it obvious that an even number of intervals 


( odd number 

of 

points ) are required 

for the 

application 

of Simpson's 

1/3 

Rule . The 

smal 1 er 

we 

choos e 

Az more 

accurate are 

the 

results . 





In 

equation (3.9) 

, surface 

is 

the 

function of 


height of the storage unit . Its values at each of the 
node points are obtained from the steady state temperature 
distribution obtained previously . Finally a computer 
programme has been prepared to evaluate above integral 


The 

integral evaluates the 

value of 

heat 

generated 

and 

this 

value is compared 

with the 

original value 

of 

( q’ 

’ ’x volume ) the heat 

generated 

in a 

storage unit 

- 

These 

values should match 

with each 

other 

. However due 

to 


numerical calculations there is a small difference between 
these values in the results obtained by computer 
programme . 

The variation of temperature with height and 
diameter, are shown in Tables 3.1 and 3.2 . Graphically these 
variation are shown in Fig. 3.3 and 3.4 . 
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TABLE 3.1 


Variation of temperature with radius at half the height 


RADIUS 


TEHPERATURE ( deg. cent. ) 



0.00 


389.40 


289.12 


255.62 


238.86 


0.01 


389 . 40 


289.12 


255.62 


238.86 


0.03 


387.28 


287.01 


253.51 


236 .75 


0.04 


383.16 


282.87 


249 .37 


232 .61 


0.05 


0.06 


0.08 


0.09 


360.25 


349.80 


276.86 


243.36 


235.64 



226.60 


218.88 


209 .67 


249 . 48 


0.10 


338.34 


238.01 


204.50 


187.73 


0.12 


0.13 


0.14 


0.15 


0.17 


0.18 


313.35 


300 . 28 


287.09 


261.04 


225.76 


212.99 


199 . 90 



110.32 
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TABLE 3.2 

Variation of temperature with height at zero radius 


HEIGHT 

m . 


TEMPERATURE 

( deg. cent. ) 



h = lOUm'^K"^ 

h= 15Um“^K"^ 

h =20Uin 

0.00 

211.37 

133.50 

108.33 

95 . 95 

0.14 

319.50 

238.51 

211.91 

198.67 

0.29 

362 .46 

273.08 

243.05 

227 .94 

0.43 

379.23 

284.12 

252.06 

235.94 

0.57 

385.73 

287 .61 

254 .64 

238.10 

0.71 

388,24 

288.71 

255.37 

238.68 

0.86 

389.16 

289 . 05 

255.58 

238.83 

1.00 

389 . 40 

289.12 

255.62 

238.86 

1 .14 

389.16 

289 . 05 

255.58 

238.83 

1.29 

388.24 

288.71 

255.37 

238.68 

1 

1.43 

385.73 

287 .61 

254.64 

238,10 

1.57 

379 . 23 

284.12 

252.06 

235.94 

1 

1.71 

362.46 

273.08 

243.05 

227.94 

1.86 

319.50 

238.51 

211.91 

198.67 


2.00 


133.50 


108.33 


95.95 










TEMPERATURE (d 



0.00 0.05 0.10 0.15 0.20 

RADIUS (meter) 

Fig. 3.3 Variation of temperature with radius for different 
heat transfer coefficients 
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Chiapt,er 4 , 

TRANSIENT ANALYSIS 

4-.1 INTRODUCTION 

The design o£ storage vault utilises the decay- 
heat o£ the waste and engineered design of the vault and 
the stack induces air £low by natural draught to maintain 
permissible temperatures . Since this type of 

storage utilises passive cooling by naturally induced 

draught to the concrete wall structure , its effectiveness 


is certain . The 

maximum credible 

accident 

which 

could 

severely 

affect 

the 

cooling ia 

one 

in 

which 

the 

stack 

collapses 

as 

itself and 

blocks 

the 

exhaust 

plenum 

, thus 

blocking 

the 

air passage . 

Such 

conditions 

could 


presumably continue till the plenum is unblocked . A 
theoretical analysis of the maximum credible accident has 
been carried out . It should be noted that this 
analysis is based on a highly pessimistic approach .It 
has been assumed that both the exit and entrance 
plenums of the vault are blocked . Further the 
contribution of the vault walls to the dissipation of 
heat to the soil has not been considered . The actual 
temperatures are expected to be much lower than 
projected . 

Provisions have also been incorporated for forced 
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draught circulation of air employing centrifugal blowers in 
case of failure of the natural draught cooling , the 
probability of which is very remote . Arrangement for 
circulation of air in the annular space between primary 
and secondary vault are also provided . 

To estimate the consequences of such accidental 

conditions transient analysis of a storage unit is carried 
out . In this the initial conditions are same as those 

reached in steady state analysis . The use of this analysis 
is to get the time after which critical temperatures are 
reached . Curves have been plotted showing the rate of 
temperature rise . 


4-.2 METHOD OF ANALYSIS 


The basic equation of transient heat transfer is 


2 2 q ' 

1 T _ S I ^ 1 ^ T ^ S I ^ ^ 

a S~T ^2 V s v ^2 K 

^ r o z 


(4.1) 


where , 


a. = Thermal diffusivity 
T = Time 
T = Temperature 
r = Radius 
2 = height 

q'’*= Rate of heat generation 

K = thermal conductivity of the storage unit material 
Again finite difference scheme will be used to 
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'ambient ~ Ambient air temperature at time t + At 
j = Temperature at node (i.j) at time t 


T. . = Temperature at node Citj) at time t +At 

» J 

The value of T^^^ , obtained from (4.2) for every 
node point is put into (4.3) Uhich gives the value of 
^ambient ' temperature of ambient air after 

time At . Because of a large number of node points 
involved in the network a computer programme is prepared . 
Uhich gives the temperature after any time T . 


4.3 RESULTS 


Above calculations are done for three values of 
h the convection heat transfer coefficient 

1. h = 10 U m"^K"^ 

2. h = 5 U m"^K~^ 

3. h = 0.0 

to get a wide range of results . Third one is the worst 
case in which there is no convection heat transfer 
between the storage unit and surrounding air . The variation 
of these temperatures with time are shown in Figures 4.1 , 


4 . 2 and 4.3 



0.00 


Variation 

and 


50.00 100.00 150.00 200.00 

TIME (minute) 

of atmospheric temperature with time 
convection heat transfer coefficient (h) 



TEMPERATURE (deg. cent.) 



Fig. 4.2 Variation of temperature with radlua 

of the cylinder and time for h = 10 U m 



MAXIMUM j TEMPERATURE (de 



TIME (hour) 


Figure 4.3 : variation of maximum temperature 

with time for h = 10 U 
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Chapter 5 


CONCLUSIONS 


1. 

It is 

necessary to 

store 

the 

solidif ied 

wast e 

for 

an interim 

period of 

time 

to 

dissipate 

the 


decay heat to avoid melt down . 

2. Constant surveillance of the waste product and 
its container should be done until the major portion of 
decay heat is dissipated . 

3. Uater is a better coolant than air during 
interim storage because of its high heat capacity ; but 
due to corrosion problem it is not used . 

4. Actinides are relatively less important within the 
time frame of interim storage , because of their long half 


lire . 


5. 

The 

storage vault is 

to be designed on 

thermal 

as well 

as 

structural consideration . 




6. 

The 

pressure loss in 

the cross 

flow 

of 

air 

through 

the 

Storage unit is not significant 

. It 

is 

only 

about 1 

percent of total head 

developed 

However 

the 

pressure 

loss 

in the stack is 

around 60 k 

giving 

a 

stack 

height 

of 100 meter . 





7. 

Out 

of three cases of 

heat transfer , 

viz . 

flow 

across 

the tube bundles , free 

convection 

heat 

trans f er 

from a 

bare . 

cylinder and free 

convection 

through 

an 

annular 

space 

, the convection 

heat transfer 

coefficient 

is 
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lowest in the third case . However the radiation heat 

transfer coefficient is around 160 % of convection heat 

transfer coefficient in case of annular space . If it is 

possible to have perfectly black surfaces , then the 

radiation heat transfer coefficient becomes 250% of 

convection heat transfer coefficient in this case 
calculations show that the heat transfer coefficient in 

case of annular space can never exceed the same in case 

of flow across tube banks . 

8- In the steady state analysis the maximum 

temperature always occurs at the centre of cylinder for 

the whole range of convection heat transfer coefficient . 

9. Transient temperature distribution shows that the 

radial temperature distribution becomes flatter with the 

passage of time and the ambient temperature in the 

storage vault goes on rising . Calculations show that the 

time after which critical temperatures are expected to 
reach may vary from 3 to 15 hours depending upon the 

convection heat transfer coefficient in the building . 
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